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Abstract 

We used a local susceptibility approach in extensive polarized neutron diffraction studies of the spin liquid Tb2Ti2 07- For a 
magnetic field applied along the [110] and [111] directions, we found that, at high temperature, all Tb moments are collinear and 
parallel to the field. With decreasing temperature, the Tb moments reorient from the field direction to their local anisotropy axes. 
For the [110] field direction, the field induced magnetic structure at 10 K is spin ice- like, but with two types of Tb moments of very 
different magnitudes. For a field along [111], the magnetic structure resembles the so-called "one in-three out" found in spin ices, 
with the difference that all Tb moments have an additional component along the [111] direction due to the magnetic field. The 
temperature evolution of the local susceptibilities clearly demonstrates a progressive change from Heisenberg to Ising behavior of 
the Tb moments when lowering the temperature, which appears to be a crystal field effect. 
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J.. Introduction 

■ Recent polarized neutron studies of the magnetic distri- 
bution in some ferromagnetic and paramagnetic materials 
with the T/13P4 cubic structure have shown that the mo- 
ment induced by a magnetic field at equivalent crystallo- 
graphic sites may be very different |l|2|3j . Such a behavior 
can arise when the local symmetry of the magnetic ion is 
lower than the overall symmetry of the crystal. The mo- 
ment induced on each magnetic ion by the internal or ex- 
ternal magnetic field depends on the orientation of the field 
with respect to the local symmetry axis. This effect can 
be described by attributing to each magnetic atom a site 
susceptibility tensor Xij which gives the magnetic (linear) 
response of the ion to the applied magnetic field [4]. The 
symmetry of the Xij tensor is the same as that of the Uij 
tensor describing the thermal motion of atoms. The com- 
ponents of the Uij tensor represent the mean square atomic 
displacement parameters (ADPs). In the linear approxima- 
tion, i.e. for relatively weak magnetic fields, one can in- 
troduce atomic suceptibility parameters (ASPs) in a way 
analogous to the ADPs. The response of an atom to a mag- 
netic field can then be conveniently visualized as a mag- 
netization ellipsoid constructed from the six independent 



ASPs in much the same way as thermal ellipsoids are con- 
structed from the ADPs. In the absence of local anisotropy 
the magnetic ellipsoids reduce to spheres with their radii 
proportional to the induced magnetization, but in many 
cases anomalous (elongated or flattened) ellipsoids can ap- 
pear. The presence of such anisotropic magnetic ellipsoids 
accounts for the anomalous magnetic behavior mentioned 
above. 

R2Ti207 pyrochlore compounds represent another type 
of a cubic structure with space group Fd 3m where the local 
symmetry of the magnetic ion is lower than the overall sym- 
metry of the crystal. In R2Ti2 07, the rare earth ions occupy 
a 16d site with a local trigonal symmetry 3m and are sub- 
ject to a strong crystal field interaction. Both the R and Ti 
atoms independently lie on a pyrochlore structure, a face- 
centered cubic lattice of corner sharing tetrahedra which 
produces phenomena of geometrical frustration. Depending 
on the nature of the magnetic rare earth ion in these materi- 
als, the ground state can exhibit long range magnetic order 
|5|6j or spin ice physics [7 8j, and in the case of Tb2Ti2 07it 
is a highly correlated quantum disordered state known as 
a spin liquid [9 10 llj. In Tb2Ti207 [12 13j the anisotropy 
is weaker than in canonical spin ices, and ferromagnetic 
(F) and antiferromagnetic (AF) first neighbor interactions 



nearly compensate. This makes Tb2Ti207 extremely sensi- 
tive to external perturbations, which leads to a rich variety 
of magnetic ground states, from spin liquid [9j to antifer- 
romagnetic order under pressure [14j and/or in magnetic 
field |15|16j . We undertook a systematic study of the field 
induced magnetic orders in the spin liquid Tb2Ti207 in a 
wide temperature (0.3 < T < 270 K) and field {0 < H < 7 
T) range, by combining polarized and unpolarized neutron 
diffraction on a single crystal. Unpolarized neutron results 
will be published elsewhere |17|18 j ; we describe here in de- 
tail our polarized neutron results based on the local suscep- 
tibility approach. We show that this approach is extremely 
eflacient in the treatment of polarized neutron data as it al- 
lows an universal description of magnetic structures exist- 
ing in large temperature (5-300K) and field (1-7 T) ranges, 
regardless of the field direction, with only 2 parameters. 

2. Experimental 

Neutron diffraction studies were performed on the 
diffractometer Super-6T2 [19j and 5C1 at the ORPHEE 
reactor of the Leon Brillouin Laboratory, CEA/CNRS 
Saclay. Polarized neutron fiipping ratios were measured on 
Super-6T2 using neutrons with = 1.4 A obtained with 
a supermirror bender (Pq = 0.98) and completed by mea- 
surements on 5C1 (An = 0.84 A) obtained with a Heus- 
sler polarizer (Pq = 0.91). The programs CHILSQ of the 
Cambridge Crystallography Subroutine Library [20j and 
MEND [21] were used for the least squares refinements on 
the fiipping ratios and maximum entropy reconstruction 
respectively. 

A single crystal of Tb2Ti207 was grown from a sintered 
rod of the same nominal composition by the fioating-zone 
technique, using a mirror furnace [17j. Prior to polarized 
neutron measurements the crystal was characterized by 
neutron diffraction at 180 K and 2 K in zero-field. A total 
of 238 refiections with sin 0/X < 0.6 A-i was measured 
at 2 K. The structure factors of 46 unique refiections were 
obtained by averaging equivalents. These were used to re- 
fine all positional parameters, the occupancy factors, the 
isotropic temperature factors and the extinction parame- 
ters. The refinement of the crystal structure was performed 
using space group Fd3m. 

3. Magnetic ellipsoids and anisotropic 
susceptibility parameters 

According to Ref. [4J atomic (local) susceptibility pa- 
rameters can be determined from polarized neutron fiip- 
ping ratio measurements if the proper magnetic symmetry 
of the crystal is taken into account. In order to investi- 
gate the temperature dependence of local susceptibilities 
of Tb2Ti207, a series of fiipping ratio measurements was 
made at 1.6, 5, 10, 20, 50, 100, 150, 200 and 270 K with a 
magnetic field of 1 T applied parallel to the [110] direction. 
(In fact there was a misalignment of about 5 degrees be- 



tween the [110] direction and the magnetic field direction. 
The exact orientation of the magnetic field with respect 
to the cubic crystal axes was taken properly into account 
in the final data analysis.) The results were interpreted in 
terms of a model which assigns a site susceptibility ten- 
sor X to each crystallographically independent site [4J . The 
point group symmetries of the atomic sites of the paramag- 
netic group F(&m were used to constrain the site suscepti- 
bility tensors. In the cubic axes, the symmetry constraints 
for a magnetic atom occupying the 16d site in the FdSm 
group imply: xii = X22 = X33 and Xi2 = Xis = X23. Thus, 
only two independent susceptibility parameters xii and Xi2 
need to be determined regardless of the field direction. 

For each temperature these two independent components 
Xii and X12 were determined using the least squares refine- 
ment program CHILSQ [20]. The refinement carried out 
with the 102 fiipping ratios measured at 270 K with the 
magnetic field = 1 T applied parallel to [110] gave Xii = 
0.056(5) /i^/T and X12 = 0.005(6)(11) /i^/T with a good- 
ness of fit x^= 1-02. For convenience, the susceptibility com- 
ponents are given in /i^/T, which allows an easy compari- 
son with the results of a conventional least squares refine- 
ment based on the localized magnetic moment model. The 
fact that the non-diagonal term X12 is equal to zero (within 
error bars) shows that the magnetic ellipsoids at 270 K re- 
duce to spheres of diameter 0.056 /i^/T. It means that the 
moments induced on all Tb atoms in the unit cell have the 
same magnitude and are collinear with the magnetic field 
regardless of the field direction. This behavior resembles 
that of a Heisenberg paramagnet at high temperature. The 
situation is quite different at low temperature where the re- 
finement carried out with the 423 fiipping ratios measured 
at 10 K in the same magnetic field gave xii = 0.939(15) 
/i^/T and Xi2 = 0.535(11) /i^/T with the goodness of 
fit x^=2.92 . One can see that, besides a considerable in- 
crease of the diagonal term xii refiecting the Curie- Weiss 
behavior, the non-diagonal elements Xi2 become very large, 
showing that the principal axes of the magnetic ellipsoid 
are different in length and do not coincide with the cubic 
crystal axes. It is easy to show that the trigonal local sym- 
metry implies that the main principal axes of the ellipsoids 
lie along the four [111] axes and their lengths are given by: 
X||=Xii+2xi2 and x±=Xii-Xi2. 

The evolution of the magnetic ellipsoids with temper- 
ature is illustrated in FiglH One can see that the shape 
and the mutual orientation of the ellipsoids refiect the cu- 
bic symmetry of the structure. The ratio of the length of 
the main axes, X|| along the local anisotropy (trigonal) axis 
[111] and x± perpendicular to it, is strongly temperature 
dependent and the ellipsoid elongation, which is due to an 
enhancement of the Xi2 component, increases markedly at 
low temperatures. Magnetization along the easy and hard 
local axes at 10 K differ by a factor of 5, which shows that 
the local anisotropy progressively evolves from Heisenberg 
to Ising character when the temperature decreases. 

Once the site susceptibility parameters are known, one 
can easily calculate the magnitude and the direction of the 
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Fig. 1. Top: Magnetic structure (arrows) and magnetic ellipsoids measured in the field H=l T || to [110] at 10 K (a), 100 K (b) and 270 K (c). 
Ellipsoids were scaled by temperature to compensate the Curie- Weiss behavior. Bottom: Two dimensional (2D) sections of the magnetization 
density component mz parallel to the field (red contours) and three dimensional (3D) isosurfaces of the magnetization density mz (blue). 
The level of the isosurface selected on each figure corresponds to 25% of the maximal density in the cell found in the reconstruction. The 
maximum values of densities are scaled by temperature. 



moment induced on each Tb atom by a magnetic field ap- 
plied in an arbitrary direction since the site susceptibility 
tensor relates the vectors M on each Tb ion and H: M = 
X H. The magnetic structures at 10, 100 and 270 K re- 
sulting from the refinement are shown in FigHJ There is a 
progressive evolution of the structure from collinear, with 
the same magnetic moment on all Tb ions, to non-collinear 
with two types of Tb ions having a high and a low moment. 
This corresponds to the separation of Tb ions into so-called 
a and /^-chains |23|24j . Moments in a-chains have their lo- 
cal < 111 > easy axis close to H (35.3°), whereas moments 
in /3-chains have their easy axis perpendicular to H. The 
non-collinear magnetic structure induced at low tempera- 
ture by the field can be seen as a direct consequence of the 
strong local anisotropy, the moments remaining close to the 
elongated axes of the ellipsoids whatever the field direction. 

To verify the validity of the local susceptibility approach, 
a second series of fiipping ratio measurements was per- 
formed with a field of 1 T parallel to the [111] direction 
at different temperatures. Refinement on 144 fiipping ra- 
tios at 10 K gave xii = 0.95(5) ijlb and Xi2 = 0.55(4) ijlb 
with x^=7.1. The Xij parameters for the [111] field direc- 
tion coincide within experimental error with those obtained 
for the field along [110]. The magnetic structure at 10 K 
resulting from the refinement is shown in FigO One can 
easily recognize the so-called "one in-three out" spin con- 



figuration described earlier in Ref. [23] , with the difference 
that all Tb moments have additional components along the 
[111] direction due to the magnetic field. These results con- 
firm that the polarized neutron data collected with differ- 
ent field orientations can be interpreted within the param- 
agnetic group FcSm using only the 2 susceptibility param- 
eters allowed by symmetry. 




(a) (b) 

Fig. 2. a) Magnetic structure and magnetic ellipsoids measured in a 
field H=l T II to [111] at 10 K. b) The corresponding 2D map and 3D 
isosurface of the magnetization density m^, see captions to Fig. 1. 



4. Magnetization densities 

Our refinements show a significant difference between 
tiie magnetic moments in a and /3-chains and it is useful to 
check to what extent this result can be substantiated by a 
model-free analysis of our data by the maximum entropy 
method (MEM). This method has been shown to give much 
more reliable results than conventional Fourier syntheses, 
by considerably reducing both noise and truncation effects 
[22] . In order to carry out the MEM reconstructions the 
magnetic structure factors of all data sets were extracted 
from flipping ratios using the SORGAM program [20j. 

As shown in Ref. [2j when the local symmetry of the mag- 
netic ion is lower than the overall symmetry of the crys- 
tal, the high symmetry paramagnetic group Fd3m cannot 
be used in the magnetization density reconstruction. For 
the field oriented along the two- fold axis [110], for exam- 
ple, the low symmetry (orthorhombic) Fdd2 group must be 
used. This is the highest symmetry subgroup of Fd3m un- 
der which the homogeneous magnetization component mz , 
induced parallel to the two-fold axis by the magnetic field, 
is invariant. An important consequence of the symmetry 
reductions is that the number of independent reflections 
to be measured increases considerably compared to that of 
the paramagnetic group. A typical number of reflections 
used in the reconstruction varies from 100 to 400, depend- 
ing on the degree of local anisotropy and on the magni- 
tude of the flipping ratios. The magnetization density dis- 
tribution was discretized into 51x51x51 sections along a, 
b and c respectively. Then the magnetization density com- 
ponent parallel to the field {mz) was reconstructed using 
a conventional uniform (flat) density [22]. Results of the 
reconstruction are presented in FigHja-c. For each temper- 
ature, the figure shows both a 2D section of the magne- 
tization density in the [111] plane and a 3D isosurface of 
the density. One can see that at 270 K the magnetization 
isosurface is nearly the same on a and /3-chains, while at 
10 K magnetization is practically absent on the /^-chains. 
Thus the reconstruction confirms the presence of a large 
difference of mz components of the Tb moment in a and 
/3-chains at low temperature, in spite of the fact that the 
maximum entropy method is strongly biased against a dis- 
equilibrium of magnetic density in the unit cell. We note 
also that sections through a Tb site at 10 K show that the 
magnetization density on a -chains is elongated in the [111] 
directions (FigUja, bottom). This might be a consequence 
of the strong crystal field anisotropy at low temperature. 
We stress that, in the above figures, we only present the 
mz component of the density as the flipping ratio measure- 
ments can only yield the density component collinear with 
the field. On the other hand the presence of a high degree 
of anisotropy in the local susceptibility implies a consider- 
able non-collinearity on the magnetization distribution at 
low temperature. To account for this, more sophisticated 
neutron experiments are needed [25] . We suggest that the 
non-collinear part of the magnetization density can proba- 
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Fig. 3. Longitudinal and transverse local susceptibilities x\\ cind x± 
versus T. Square and circle symbols correspond to a field along [110] 
and [111] respectively. Filled triangles show the bulk susceptibility 
derived from the neutron data, x^^^=^/^ X\\ +2/3 x±^ open tri- 
angles the bulk susceptibility from Ref. [13 . Lines are calculations 
using the CF parameters of Tb2Ti207 (see text). 

bly be extracted by combining flipping ratio and integrated 
intensity measurement in a magnetic field, but this needs 
to be verified experimentally. 

A similar maximum entropy reconstruction (using a trig- 
onal symmetry subgroup of Fd^m) was performed on flip- 
ping ratios measured with the field parallel to [111]. The 
resulting magnetization density component mz obtained 
at lOK in 1 Tesla is shown in Fig|2l It gives an evidence 
that Tb ions having their local axis parallel to the field 
have a very high magnetic moment, while the rest of the 
Tb ions have very small moments, in agreement with the 
biased "one in-three out" structure obtained by the refine- 
ment (bottom FigUja). 

5. Longitudinal and transverse local 
susceptibilities 

Figure[3]shows that the ratio of the main axes of the ellip- 
soids X||/x± increases with decreasing temperature. This 
reflects the increasing CF anisotropy as temperature de- 
creases. Indeed, at high temperature, many CF states are 
populated and the anisotropy is close to that of the free ion 
(weak or Heisenberg-like), whereas at low temperature, the 
lowest doublets alone are populated, which have a strong 
Ising behaviour. The local susceptibilities obtained with 
the field along [110] and [111] are in a very good agreement, 
which confirms that this approach is valid regardless the 
field orientation. 

For comparison with the experimental data, we used the 
CF parameters of Tb2Ti207 derived from the inelastic neu- 
tron spectra [13j. The parameters account properly for the 
trigonal symmetry of the Tb environment and allow to 
compute the magnetization induced by a magnetic field ap- 
plied parallel or perpendicular to the local < 111 > axis. 
Interatomic exchange and dipolar interactions were taken 
into account in the molecular field approximation via a mi- 
croscopic constant A. The self-consistent calculation, valid 



down to 10 K, yields a constant A = -0.35 T/ jiB^ close to 
the value derived in Ref. [13j. The result of the calculations 
agrees well with the experiment, see Figl3l 

In the local susceptibility approach, the atomic site sus- 
ceptibility tensor accounts only for the linear paramagnetic 
response of the moments to an applied field. The range of 
fields and temperatures satisfying the condition of linear re- 
sponse can be estimated from the bulk magnetization data 
[13j: the Xij parameters obtained in our refinement allow 
the general behavior of Tb2Ti207 to be described above 10 
K in fields up to 7 T applied in an arbitrary direction. 

6. Conclusion 

We determined the field induced magnetic structures in 
Tb2Ti2 07 by polarized neutron diffraction, in a wide range 
of temperatures and fields, applied in the [110] or [111] di- 
rections. The [110] low temperature ground state shows a 
non-collinear magnetic structure that resembles the local 
spin ice ground configuration, but with Tb moments of very 
different magnitudes. The [111] low temperature ground 
state is close to the structure "one in-three out" found in 
spin ices, with the difference that all Tb moments have an 
additional component along the [111] direction due to the 
magnetic field. This is attributed to the finite CF anisotropy 
of Tb2Ti207 which is much weaker than in model spin ices. 
Polarized neutron diffraction and the local susceptibility 
approach provide an easy way for studying field induced 
structures in pyrochlore compounds. They provide a uni- 
versal description of the magnetic structures in Tb2Ti207 
in a large range of fields and temperatures, with only two 
local susceptibility parameters regardless of the field direc- 
tion. 
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